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Recently, interesting data concerning dN^/drj in Au-Au collisions [rj = — lntan((9/2)] 
with centrality cuts have been reported from the PHOBOS Collaboration. In most treatment 
these data are divided by the number of participants (nucleons) in collisions. Instead of this 
£f~) • method, we use the total multiplicity N c h = J {dN c ^/ drj)drj and find that there is scaling 

phenomenon among (N c i 1 )~ 1 dN c i 1 /drj = dn/drj with different centrality cuts at y / ijvjv = 130 
GeV. To explain this scaling behavior of dn/drj, we employ a stochastic approach using 
the Ornstein-Uhlenbeck process with two sources. A Langevin equation is adopted for this 
explanation. Moreover, comparisons oidn/drj at ^/snn = 130 GeV with that at -^/sjvjv = 200 
GeV are made, and no significant difference is found. A possible method fot the detection 
of the quark-gluon plasma (QGP) through diVch/dr? is presented. 

O 

§1. Introduction 



Recently, interesting data from the PHOBOS Collaboration on dN^/dn [rj 



— lntan(0/2)]*)in Au-Au collisions at ^snn = 130 GeV have been published. ^ The 



authors of Ref. 1) calculated the quantity 



(CT^f = /(<Ari " ) ' Ar ~ l "" ) ' (M) 

where N part and iV co ii represent the number of participants (nucleons) and the num- 
ber of collision particles in Au-Au collisions. The quantity in Eq. (IT) depends on 
the centrality cuts. The intercept /(iV part , N co n, n = 0) is an increasing function of 

(AW)- 



Here, 



yT + m 2 /pt + sinh 2 rj + sinh rj 
2 E — p z 2 yl + m 2 /pi + sinh 2 rj — sinh rj 



1 m^ = Im 



= tanh -1 H-j^J » - Intan(0/2) = rj . 



1 P + Pz dn p dn p cosh rj 

rj — — In and — = — , where — = — ^^^^^^=^^^^^^= 

2 p-p z drj E dy E y / l + m 2 / p 2 + sinh 2 ~ 
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In this paper, instead of Eq. (1-1), we consider the physical quantity 

1 (i/Vf-h dn 
N ch dn drj ' 



(1-2) 



where = J '(dN c ^/ 'drj) drj and J(dn/dn)dn = 1. In Fig. 1, three sets of values for 
dn/dr] are shown. They suggest that there is scaling behavior among the different 
sets of dn/dn with different centrality cuts. Thus, dn/dn can be considered a kind 
of probability density, because the variable n is a continuous variable. This fact 
probably implies that a stochastic approach is appropriate in analyses of dn/drj.*^ 

The content of the present paper is as follows. In §2, we examine the dn/dr] 
scaling. In §3, a stochastic approach is considered as one possible explanation of 
dn/drj scaling. In §4, concrete analyses employing Gaussian distributions obtained 
from the Ornstein-Uhlenbeck (O-U) process are presented. In the final section, 
concluding remarks are given. In the Appendix, the Fokker-Planck equation for the 
O-U process is considered. 

§2. Confirmation of dn/dr] scaling 

It is worthwhile to confirm whether dn/dr] scaling holds. Using the intercepts of 
dn/dn at n = in Fig. 1, 



dn 
drj 



= c^0.129±0.005, (2-1) 

T]=0 



we can obtain a relation between N c h and (iV part ) as follows. The intercept at rj = 
can be parameterized as 



1 dN ch 



0.5(iV part ) dn 



= A(N part ) a , (2-2) 

T]=0 



*^ From studies of multiparticle dynamics in high energy physics, we have learned that the 
probability distributions P(n, (n)) are functions of n and (n). It is known that the KNO scaling 
functions 2 ^ 

lim (n)P(n) = ip(z = n/(n» , 

are described by solutions of various Fokker-Planck equations. 3 - ) ' 4 ^ Moreover, this stochastic descrip- 
tion is also seen in QCD. For example, Dokshitzer has calculated the generalized gamma distribution 
in QCD 5) 

P(z = n/(n}) « ^ ^ z > exp [-(£>*)"] , 

where z is the KNO scaling variable, fj, is given by 1 — 7 = 1/fj,, D is a parameter, and 7 is the 
anomalous dimension in QCD. This is a stationary solution of the following Fokker-Planck equation 



ill' i) ■ 1 ■•»- 

~dt ~ ~~d~z 
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Fig. 1. (a) and (b) Two sets of dnjdr\ with different centrality cuts. 1 ' Data with 15-25% are plotted 
in both figures for the sake of comparison, (c) To examine dn/d-q in the fragmentation region, 
a log-linear plot is given. 



where A = 2.16, and a = 0.064. Thus we obtain the relation 



0.5(iV part ) 1 dN ch 



N ch 0.5(iY part ) drj 



O-^Part) A{NpaTt) a . (3.3) 



„=0 *<* 



Equation (2-3) is examined in Table I and Fig. 2. That Eq. (2-3) holds approximately 
among 6 centrality cuts reflects the dn/drj scaling, in particular, in the central region. 

§3. A possible explanation of dn/drj obtained using a stochastic 

approach 



It is well known that the rapidity (y rj) is a kind of velocity. Moreover, there 
are leading particles in the beam and in target nuclei that collide with each other. 
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Table I. Empirical examination of Eq. (2-3). 



centrality (%) 


2c x N ch 


A x <iV part > 1+a 


35-45 


2 x 0.129 x 1056 


T — TTTTTrn 

2.16 x 93 1 + U UM 




272.6 ± 14.6 


268.5 ± 17.7 


25-35 


2 x 0.129 x 1582 


2.16 x 135 i+0 uti4 




408.1 ± 21.9 


399.1 ± 28.4 


15-25 


2 x 0.129 x 2270 


2.16 x 197 1+U ' UM 




585.7 ±31.4 


596.7 ±45.8 


6-15 


2 x 0.129 x 3199 


2.16 x 270 i+U UM 




825.2 ±44.6 


834.5 ± 67.9 


0-6 


2 x 0.129 x 4070 


2.16 x 340 i+u ' UM 




1050 ± 57 


1066 ± 90 




For example, nucleons in gold atoms in RHIC experiments collide with each other, 
and thereby lose energy and emit various particles. Since we have to treat large 
numbers of particles (e.g., lk-lOk), a stochastic approach seems to be appropriate 
and to offer a simpler description than Monte Carlo approaches.*) 

To describe dn/dr] scaling with the leading particle effect and fluctuations in 
rapidity space, we assume the following Langevin equation 8 ) " 10 ) for the rapidity 
variable:**^ 

^ = -7V + /w(*). (3-1) 

Here t, 7 and f w (t) are the evolution parameter,*) the frictional coefficient and a 
white noise term, respectively. In our treatment, we assume that particles are 

*^ It should be noted that the transport approach is another useful method (see Refs. 6) and 

?))• 

**) In classical mechanics, Eq. (3-1) corresponds to the equation 

m— = -mfv + m/ w (t) , 

where v is the velocity. 

*' As an alternative interpretation, t may be related to the number of collisions among the wee 
partons and produced particles. 
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produced at ±y max at t = 0. This picture takes into account leading particle effects. 
Using the assumption y(0) = ±y max , we obtain the solution 

y(t) = ±y max e- 7 * + e" 7 ' /* e^ s f w (s)ds . (3-2) 

Jo 

The average and variance of y{t) are calculated as 

E[y(t)} = ±y max e- 7 * , (3-3) 



E[(E[y(t)\ - y) 2 ] = °- (l - e" 2 *) , 
where we use the expression 

(U(t)U(s))=^5(t-s) 



(3-4) 



(3-5) 



for the white noise, where a 2 is the variance. It is known that the distribution 
function for y(t) is given by a Gaussian distribution with the above average and 
variance. The probability density with V 2 (t) = (o 2 /2j)(l — e~ 2lt ) is given by 



P(y, 2/max, t) 



1 



cxp 



+ exp 



(V + V* 



2V 2 (t) 

(y - ymaxe" 7 *) 2 



2V 2 {t) 



(3-6) 



The connection between Eq. (3T) and the Fokker-Planck equation for the O-U pro- 
cess is given in the Appendix. 

In Fig. 3(a), we depict a simplified picture of heavy-ion collision. Our assump- 
tions for the leading particle effect are equivalent to the assumption P(y, y m ax, t = 
0) = 0.5[<5(y — y m ax) + £(y + ymax)]- I n other words, in this model we make the simple 
assumption that there are two sources of particles £tt t — ; located at iymax 

and 

producing 0.5-/V c h particles each. The evolution of P(y, y m ax 5 i) given in Eq. (3-6) is 
shown in Fig. 3(b). 



§4. Analyses of dn/dr] by means of Eq. (3-6) 

Making use of Eq. (3-6), we can analyze dn/dr] shown in Fig. 1. In our calcula- 
tion, as most produced particles are not specified, we assume that y « r\ in Eq. (3-6). 
Our results are shown in Fig. 4 and Table II. In Fig. 5, we examine whether or not 
the variance V 2 (t) and the quantity p = 1 — e -27 * depend on the centrality cuts. As 
is seen in Figs. 3 and 4, the scaling behavior among the sets of dn/dn at ^/snn = 130 
GeV is explained by Eq. (3-6), with small changes in the variance V 2 (t). The values 
of V 2 (t) depend on the distribution in the fragmentation region [— n max < n < —4 
and 4 < n < rf max \. It can be said that the scaling behavior is explained fairly well 
by the O-U process with two sources, one at the beam (ye or y max ) rapidity and one 
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target 
fragmentation 
region 



f«0 



(b) 



f«0 




t — > finite 




Fig. 3. (a) Simplified picture for A-A collisions. The thin lines represent partons, and the black 
circles represent nucleons. (b) Evolution of P(y, 2/ ma x, t) in Eq. (3-6) with two sources at t/ max 
and -j/ max . 



at the target (yx or — y max ) rapidity. We have confirmed that an O-U process with 
a single source is not capable of explaining the scaling behavior. *)>**) 

Table II. Parameter values obtained in our analyses using Eq. (3-6) with two sources. The evolution 
of P(y, 2/max, t) in Eq. (3-6) is stopped at minimum values of \ 2 - Here, Sp = 0.006 0.004 

. (n.d.f means the number of degree of freedom.) 



Fig. 4 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


centrality (%) 


45-55 


35-45 


25-35 


15-25 


6-15 


0-6 


V 


0.872±(5p 


0.875±(5p 


0.878±5p 


0.882±5p 


0.886±(5p 


0.888±(5p 


V\i) 


3.83±0.27 


3.61±0.21 


3.23±0.16 


3.00±0.13 


2.72±0.10 


2.47±0.08 


(AWt) 




93 


135 


197 


270 


340 


A ch 


662±10 


1056±16 


1582±23 


2270±34 


3199±49 


4070±63 


c 


0.125 


0.127 


0.128 


0.130 


0.132 


0.131 


X 2 /n.d.f. 


8.61/51 


7.63/51 


5.88/51 


5.35/51 


3.57/51 


3.82/51 



and 



V2TV2(i) 



exp 



The intercepts for the sets of dn/drj are calculated using the expression 



exp 



(ymaxe 7f ) 2 

2V 2 {t) 



(4-1) 



The results are listed in Table II. They are almost the same as the values in Fig. 1. 

Here we should carefully examine the values of V(t) in Table II. The slight 
change reflects the discrepancies in the fragmentation region. As seen in Fig. 1(c), 
there are small differences in the sets of dn/drj for |r/|^4 between the 0-6% centrality 
cut and the other centrality cuts. To explore the differences more carefully, we need 



^ For an explanation with the single source at yo w 0, we can use Eq. (A-2) in the Appendix. 
For the centrality cut 0-6%, we obtain \ 2 = 25.76, with m/pt = 1.3 ± 0.1, which is necessary for 
the single source model. If m/pt is ignored, we obtain a worse value of y 2 • Thus we disregard this 
model. 

**' We have investigated whether or not the dip structures at r\ w can be explained by the 
Jacobian p/E, and obtained the worse values of \ 2 than those listed in Table II. This fact is probably 
related to the masses of produced particles that are not measured. 
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-6 -4 -2 2 4 6 



Fig. 4. Analyses of dn/dn using Eq. (3-6). (See Table II.) 



sets of dn/dri with smaller centrality cuts, for example 0-3% - 0-5%. *)' 

§5. Concluding remarks 



First, it can be said that there is scaling among the different sets of dn/dr] with 
various centrality cuts at ^/s^n = 130 GeV, as seen from the nearly constant values 
of c and the behavior in Figs. l(a)-(c). 



A calculation based on QCD is given in Ref. 11) as 



diV cl 



(Npart) dr) 



rj=0 



I so 



lot 



(^Q i 



Mi;) 



where a w 0.82, /Iqcd = 0.2 GeV and A = 0.25, and the centrality dependence of the saturation 
scale is Qqs- I n tne fragmentation region, this expression needs cutoff factors. 

Very recently, the results for dn/dr/ at ^/sjvjv = 130 GeV and 200 GeV with centrality 
cuts 0-5% and others obtained by the BRAHMS Collaboration 12 ^' 13 - 1 were reported. We have 
analyzed them using Eq. (3-6) and obtained an almost constant value of V 2 (t), because data in the 
fragmentation region are lacking for |?7|>4.5. 
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45-55 35-45 25-35 15-25 6-15 0-6 
centrality cuts (%) 



Fig. 5. The values V 2 (t) and p of Fig. 4 and Table II. 



Second, the scaling behavior of dnjdr\ is described by the solution given in 
Eq. (3-6) of the Langevin equation with two sources.*) (See Fig. 4 and the values of 
c in Table II.) 

Third, we can add the following fact. Very recently, data for dN^/dn with 
centrality cut 0-6% at y / s]v r /v = 200 GeV from the PHOBOS Collaboration were 
reported. 16 ) They are compared with data for dn/drj at t/snn = 130 GeV in Fig. 6. 
From this comparison, it is obvious that the scaling of dn/dr] holds between ^/snn = 
130 GeV and 200 GeV.**) 

Moreover, we can consider the dn/drj scaling from a different point of view, i.e., 
regarding the scaling property of Gaussian distributions. Using r/ rms = y/ (n 2 ) = 
\/J2 ri 2 dn/dn, we can compute the following quantity with z r = rj/r) rms : 

Irms^ = ^ = f(Zr = V/Vrms) • (54) 

In Fig. 7, we display our result for Eq. (5-1) and give ?7 rms in Table III. This scaling 
of f(z r = rj/r] Tms ) reflects the fact that the sets of dn/dn are described by a Gaussian 
distribution. In other words, the description of dn/dn using the O-U process is 
appropriate. 

This dn/dn scaling suggests that dn/dr] with the centrality cut 0-6% does not 
exhibit singular or particular phenomena related to signatures of the quark-gluon 

*-* To estimate the "thermalization time" of the QGP, Hwa has considered the Fokker-Planck 
equation for the motion of quarks and gluons in nuclei. 14 ' (See also Ref. 15), in which the Wiener 
process is considered for the problem of the thermalization of quarks and gluons.) 

Using Bjorken's picture 17 ' for the calculation of the energy density for \Arj\ < 0.5 with a 
geometrical picture of gold (R T w 6-7 fm, cro w 1-2 fm, V w ttR^(cto) w 300 fm 3 ), we obtain the 
following values: 



3 1 dN ch 



2V dr) 

1.2 GcV/fm 3 (200 GeV) 



1 GeV/fm 3 (130 GeV) , 



\A V \<0.5 

3 
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An 
dr| 



0.16 



0.12 



0.08 



0.04 - 




Fig. 6. Comparisons of dn/dr) with the centrality cut 0-6% at ^/snn = 130 GeV and 200 GeV. 
The solid curve is obtained for latter energy. Here p = 0.879 ± 0.007, V 2 (t) = 2.67 ± 0.24, 
X 2 /ri-d.f. = 0.63/22 and c = 0.126. 



■3 
1 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 






45-55% 
35-45% 
25-35% 
15-25% 
6-15% 
0-6% 



-2-10 1 2 

Zr = ^rms 

Fig. 7. Normalized distribution of dn/dz r with z r = rj/r) rrns scaling. The solid curve is obtained 



by|2-= ; 1 = (exp 
0.502, x 2 / n - d - f - = 92.9/321. 



+ exp 



=£]}, 



2.23, v; 2 (t) 



plasma (QGP).*)Of course, care must be taken when handling averaged quantities 
in statistics. Thus, while at present, it appears that the QGP is not created, it is pos- 
sible that the QGP is created but its signature is washed out by strong interactions 
between hadrons.**) 



From studies of the HBT effect at RHIC, the authors of Ref. 18) have concluded that the 
expected large radius from the QGP is not observed. 

Through the measurement of vi (flow of nucleons), some physicists have conjectured that 
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Table III. Values of r/ lrnB at ^/snn = 130 GeV 



Fig. 4 


(a) 


(b) 


(c) 


(d) 


(e) (f) 


centrality (%) 


45-55 


35-45 


25-35 


15-25 


6-15 0-6 




2.52 


2.49 


2.45 


2.41 


2.37 2.33 



To investigate particular phenomena, like the turbulence and/or deflagration in 
dN c h/dr], we need to analyze a single event with smaller centrality cut than 0-6%. 

Furthermore, event-by-event analyses using intermittency 21 ) " 26 -* and wavelets 27 ) 
are necessary to investigate the detection of QGP and a disoriented chiral conden- 
sate (DCC). For the latter case, the ratio of neutral pions ((vr )) to charged pi- 
ons ((7r ch )), (7r°)/(7r ch ), should be measured. These methods should be applied to 
dN c h/dri with smaller centrality cuts and larger particles. They seem to be capable of 
extracting useful information on QGPs and DCCs from the analysis of single events. 

Finally, we should mention the results of recent analyses of data concerning 
dN c h/dri carried out by the NA50 Collaboration. 28 ) It has been reported that there 
is 77-scaling of dNch/drj in the data obtained from Pb + Pb collisions at SPS. 28 ) This 
behavior can be described by a single Gaussian distribution. (See Eq. (A-2).) This 
fact suggests that the stochastic approach can be used to describe the behavior of 
dNffo/dr] at SPS. The intercepts at 77 max > (N milx )~ 1 dN c i l /dri = dn/dn are in range 
0.246-0.266 for 6 centrality cuts. It should be noticed that there is no dip structure. 
Their observation supports the validity of our stochastic approach. 

Note added in proof: After completion of this paper, we are informed that the rela- 
tivistic diffusion model has been used for analysis of the proton distributions in heavy 
ion collisions by Wolschin. 30 ) A similar Fokker-Planck equation is used therein. 
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Appendix A 

Solution of the Fokker-Planck Equation for the O- U Process 



The Fokker-Planck equation for the O-U process connected with Eq. (3-1) is 
given by 



dP(y, t) 



d la 2 d 2 



dy 27 dy 2 



P(y, t) . (A-l) 



the increasing value of V2 reflects the effect from the compression of nuclear matter 19 '. (See also 
Ref. 20).) 
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The solution of Eq. (A-l) with P(y, 0) = 5{y — uq) is obtained as 



P(y, t) 



y/2irV*{t) 



cxp 



(y - yoe~^ f 

2V 2 (t) 



(A-2) 



where we have used a resolution method for partial differential equations, includ- 
ing the equation of the characteristic. 29 ^ Similarly, the solution of Eq. (A-l) with 
P(y, 0) = 0.5 [S(y + y max ) + 8{y - y max )] is obtained as 



P(y, ymax, t) 



cxp 



+ cxp 



(V + Vn 



2V 2 (t) 



-jt\2 



2V 2 (t) 



(A-3) 



10 

11 

12 

13 
14 
15 

16 
17 
18 
19 
20 
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22 
23 
24 

25 

26 
27 



28) 
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Addenda 

Scaling Behavior of (N^dN^J dr] at ^I^n = 130 GeV by the PHOBOS 
Collaboration and Its Implication 

- A Possible Explanation Employing the Ornstein-Uhlenbeck Process - 

Minoru Biyajima, Masaru Ide, Takuya Mizoguchi and Naomichi Suzuki 

Prog. Theor. Phys. 108 (2002), 559-569. 

Very recently PHOBOS Collaboration has reported new analysis of dN^/di] at 
V^]\w = 130 GeV. V Thus we should add several new results based on the stochastic 
theory using their new data. 

Before our new analyses, the following re-calculations are necessary in Fig. 2 and 
Table I, because an empirical value at (-/Vp art ) = 93 and the magnitudes of error bar 
should be corrected. The method of the linear regression is used, because values of 
(-/Vpart) are computed by the Glauber model. There is no change in our concluding 
remarks, since Eq. (2.1) is explained by Eq. (2.3). 
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Fig. 2. Parameters A and a are determined by Fig. 8. Same as Fig. 2, this calculation is done 
the method of linear-regression. The cor- by new data. 1 ' The correlation coefficient 

relation coefficient (c.c.) is 0.935. (c.c.) is 0.993. 



Table I. Empirical examination of Eq.(2.3). 



centrality (%) 


35-45 


25-35 


f5-25 


6-f5 


0-6 


(iVpart) 


93 


f35 


f97 


270 


340 


N ch 


f030± ffO 


f 550 ± f 60 


3f60 ± 330 


2230 ± 230 


4030 ± 430 


c (Eq. (2.3)) 


0.f27±0.0f4 


0.f27±0.0f3 


0.f33±0.0f4 


0.f32±0.0f 


0.f34±0.0f4 



Next we have to examine Eqs. (2.1)-(2.3), using new data in Ref. 1). The 
intercept of empirical values with 6 centrality cuts is dn/d??|^ = o = c ~ 0.132 ± 0.009. 
For Eq. (2.2), we obtain the parameters A = 1.79 and a = 0.103. Eq. (2.3) becomes 
as follows 

c = °- 5( ^ part) • 1.79(AT part )°- 103 . 
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See Fig. 8 and Table IV. Because of larger c.c, we can conclude that Eq. (2.3) holds 
for new data of Ref.l), i.e., c = 0.132 ± 0.005. 

Third, we apply Eq. (3.6) to dn/dr] distribution and obtain estimated values 
given in Table V. For the z r = r)/-q rias scaling, Eq. (5.1) and the following 



Table IV. Empirical examination of Eq.(2.3) in new data. ^ 



centrality (%) 


35-45 


25-35 


15-25 


6-15 


0-6 


(iVpart) 


90 ±5 


135 ±6.0 


196.5 ±7.5 


274.5 ± 8.5 


342.5 ± 11.0 


N ch 


1000 ± 50 


1540 ± 80 


2270 ± 110 


3230 ± 160 


4100 ±210 


c (Eq. (2.3)) 


0.128 ±0.010 


0.130 ±0.009 


0.133 ±0.009 


0.136 ±0.008 


0.136 ±0.008 



Table V. New results from analysis of dN c h/drj at y / ijvjv = 130 Gev in Ref. 1) by means of Eq. (3.6), 
c(Th) = (l/y/2irV 2 (t)J exp [-(^axe" 7 ') 2 /(2V 2 (t))] and j] max = 4.9. Here, Sp = 0.004-0.006 
and 8c = 0.005-0.006. (n.d.f. means the number of degree of freedom.) 



centrality (%) 


45-55 


35-45 


25-35 


15-25 


6-15 


0-6 


?]rms 


2.51 


2.48 


2.44 


2.40 


2.36 


2.31 


p 


0.839 ± Sp 


0.845 ± Sp 


0.850 ± Sp 


0.857 ± Sp 


0.862 ± Sp 


0.866 ± Sp 


V\t) 


3.46 ±0.25 


3.33 ±0.23 


3.07 ±0.18 


2.87 ±0.15 


2.66 ±0.12 


2.46 ±0.11 


c(Th) 
X 2 /n.d.f. 


0.123 ±8c t 


0.125 ±Sc t 


0.127 ±Sc t 


0.130 ± Set 


0.131 ±Sa 


0.132 ±Sa 


2.3/51 


2.3/51 


2.2/51 


1.6/51 


2.0/51 


2.7/51 



Eq. (B-l) are applied to data of Ref. 1). The distribution by z r scaling variable is 
shown in Fig. 9. Theoretical formula is given as 



dn 1 
~dz~ r ~ v/8vrV; 2 (t) 



cxp 



-jt\2 



cxp 



-yt\2 



(B-l) 



with z max = r?max/^rms and V r 2 (t) = V 2 (t)/i]2 ms . The denominator of 
averaged rj rms among 6 centrality cuts. Figure 9 shows better z r scaling than the 
distribution of Fig. 7. It can be stressed that this z r scaling seems to be useful to 
compare dn/dn's at different energies. 2 ) 
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